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■ ABSTRACT 

O ' In the magnetically torqued disk (MTD) model for hot star disks, as proposed 

and formulated by Cassinelli et al. (2002), stellar wind mass loss was taken to be 



O 



'nI" . uniform over the stellar surface. Here account is taken of the fact that as stellar 

spin rate So (= -s/VtlR^ /GM) is increased, and the stellar equator is gravity 
p darkened, the equatorial mass flux and terminal speed are reduced, compared 

Q . to the poles, for a given total M. As a result, the distribution of equatorial 

disk density, determined by the impact of north and southbound flows, is shifted 
■ further out from the star. This results, for high So (> 0.5), in a fall in the 

disk mass and emission measure, and hence in the observed emission line EW, 
^ . scattering polarization, and IR emission. Consequently, contrary to expectations, 

^ . critical rotation 6*0 ^ 1 is not the optimum for creation of hot star disks which, in 

terms of EM for example, is found to occur in a broad peak around So ~ 0.5 — 0.6 
depending slightly on the wind velocity law. 

The relationship of this analytic quasi- steady parametric MTD model to other 
work on magnetically guided winds is discussed. In particular the failures of the 
MTD model for Be-star disks alleged by Owocki & ud-Doula (2003) are shown to 
revolve largely around open observational tests, rather in the basic MTD physics, 
and around their use of insufficiently strong fields. 

Subject headings: stars: emission-line. Be - stars: magnetic fields - stars: mass 
loss - stars: rotation - stars: winds, outfiows - polarization. 
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1. Introduction 

Be stars arc defined as ' non-sup ergiant B-type stars whose spectra have, or had at one 
time, one or more Balmer fines in emission' (CoUins 1987). Tfie pioneering researcfi on Be 
stars by Struve proposed a rotational model with emission lines from equatorial disks (Struve 
1931), but Be star disks remain enigmatic despite many decades of detailed observations 
and research (Jaschek & Groth 1982; Slettebak 1982; UnderhiU & Doazan 1982; Slettebak 
& Snow 1987; Slettebak 1988; Smith et al. 2000; Porter & Rivinius 2003). The main physics 
problems they pose are how the material in them (a) is delivered from the star; (b) becomes 
so dense; (c) acquires such high angular momentum. The answer to (a) undoubtedly lies 
in stellar radiation pressure. The first qualitative answer to (b) was the Wind Compressed 
Disk (WCD) model of Bjorkman & Cassinelli (1993). In this, the angular momentum of 
rotating wind fiow returns the matter to the equator where north and south streams collide 
and create a shock compressed disk. There are several snags with this model. Firstly, it 
does not produce high enough densities. Secondly, the disk formed has mainly radial flows 
rather than the quasi-Keplerian azimuthal flows observed. Thirdly, non-radial line-driving 
forces (Owocki et al. 1996) may cause polar rather than equatorial flow to dominate. 

A phenomenological solution to these problems was proposed and quantifled parametri- 
cally in the Magnetically Torqued Disk (MTD) model of CassineUi et al. (2002). This invokes 
a dipole-like field which steers the wind flow toward the equator and torques up its angular 
momentum on the way. The field torques up the wind fiow to Keplerian speeds or higher and 
confines the radial fiow redirecting it to be poloidal and creating a shock compressed equato- 
rial disk. The isothermal disk grows in thickness (but not in density) over comparatively long 
timescales years) which are roughly consistent with long time-scale variability of some 
Be stars (Doazan 1982; Dachs 1987; Okazaki 1997; Telting 2000), allowing a quasi-steady 
treatment. There is growing observational evidence of reasonably strong fields (hundreds of 
Gauss) in hot stars — e.g., u; Orionis with B ~ 530 ± 230 G (Neiner et al. 2003), P Cephei 
with B ~ 360±30 G (Donati et al. 2001), 6^ Orionis C with B ~ llOOilOO G (Donati et al. 
2002), though some of these are very obhque and/or slow rotating and the MTD model not 
directly applicable in its basic form. 

In the MTD treatment the stellar wind mass fiux and wind speed were taken to be 
uniform over the stellar surface. For the case of a rotating star, essential to creating a disk, 
this assumption is invalid. Rotation results in equatorial gravity darkening which reduces 
the wind mass fiux and speed there, as described by Owocki et al. (1998). In this paper, 
we evaluate the effects of this on the MTD model. We do so by generalising the basic 
quasi-steady parametric approach of MTD, but discussing in Section 5 issues concerning 
the properties of that description in relation to other theoretical and observational work on 
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the problem. The MTD paper was really the first to model the combined effects of field 
and rotation but several earlier papers had discussed disk formation by magnetic channeling 
(Matt et al. 2000), while work subsequent to MTD has variously challenged it (Owocki & 
ud-Doula 2003), and supported it (Maheswaran 2003). Insofar as there remains a degree 
of disagreement in the literature over whether the model works for Be stars, our detailed 
results should be treated with some caution, though the general trends of the effect of gravity 
darkening should be sound. 



2. Effect of Gravity Darkening on Disk Density 

Owocki et al. (1998) showed that the mass flux from a stellar surface satisfies 

Fmo{g) oc g, (1) 

while the terminal speed satisfies 

v^ig)cxg'/^ (2) 

where g is the local effective gravity. 

On a rotating sphere with stellar radius it! and angular velocity Qq, at colatitude 9, the 
nett gravity is 

GM jn^Rsiney . 

^go(l-Slsm'e), (3) 

where = ^/^J^JGM as in MTD. Strictly speaking, we should also include the effect of 
continuum scattering radiation pressure at least, which results in go = ^^(1 — Trarf), where 
Trad = ratio of the stellar luminosity to the Eddington luminosity (Maeder & 

Meynet 2000). 

It follows that the mass flux at 9 becomes 

F^oie)^Kil-Slsin^e), (4) 
where is a constant, and the terminal speed for matter from 6 is 

vM-Vooo{i-S!sm'9Y/\ (5) 
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where Vooo is the value oi Voo sit 9 — 0. 

Here we will assume the wind velocity obeys v^(r) = 1^00(1 ~ 7)'^) but with 

vUr,e)^v^{e){l--f, (6) 

r 

where (3 is assumed not to depend on 6. This basically requires that wind acceleration occurs 
quite near the star and that the field lines are roughly radial there. We want to express K 
in terms of the total mass loss rate M 



M= Fmo{0) 27r R sin eRde 
Jo 

AirR^K / {1-Sl sin^ 9) sin 0d0 
Jo 



002 

4nR'K{l-^). (7) 



So by equation (4), 

To relate this mass flux at 9 on the stellar surface to that normal to the equatorial plane at 
distance x — r/R,we follow Cassinelli et al. (2002) in parametrizing the decline of magnetic 
held, B, with equatorial plane distance according to 

B{x) = B,x-\ (9) 

where Bo is taken as uniform over, and normal to, the stellar surface, and 6 is a constant 
with a dipole field 6 = 3. Flux conservation then requires that a cross-scctional area dA of a 
flux tube arriving at the equator is related to its area dA^ at the star by BdA — BgdA^ or 

dA{x)=dAoxK (10) 

However, dA{x) = 2-KR^xdx and dAo = 2t{R^ sin 6d0 so that the relation between and x is 
given by integrating smOdO = x~'^~^^dx, to yield 

j.-b+2 

cos9 = C-j-^. (11) 
Requiring that x — > 00 for ^ — > implies C = 1 so that 

cos^(x) = l--^, (12) 
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For the particular case 6 = 3 (the dipole case) mainly discussed by Cassinelli et al. (2002), 
and which we focus on henceforth, this yields 

cos^(a;) = 1 - -. (13) 

X 

We can now obtain the mass flux F^{x) at x near the equatorial plane by using F^{x) — 
Fjno{d)dAo{9) / dA{x) and using equations (8), (10) and (13), namely 

while by equations (5), (6) and (13), the wind speed there is 

= ^ooo(l - -)^{1 - 5o[l - (1 - (15) 

and the arriving wind ram pressure Prami^) — Fm{x)vyj is 

The high density Pd{x) of cool, shock - compressed, disk (WCD) in the equatorial plane is 
then given as in the WCD model (and in MTD) by the isothermal disk (sound speed Cg) 
pressure balance expression PdC^ = Pram = Pd or by equation (16) 

Pd{So, x) = poX (1 - -f ^ _ ^gy^ ■ (17) 

Here 

which is equivalent to poc in the MTD model. Based on equation (17), we can write the 
disk density allowing for rotational gravity darkening, compared to that neglecting it (e.g., 
MTD) as 

Pd{0,x) 1-25^/3 

with the property ^{So, 1) = (1 - S^f/^/{l - 2S^/3) which is < 1 for all So, and tends to 
as — > 1. This is because the equatorial wind flow falls with increasing So, reducing 
the inner disk compression. On the other hand, if we (formally — see comments below in 
Section 3) apply equation (19) as x — > oo we would get for the disk behavior 

*(5„.oo)^^-i^. (20) 
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This is always > 1 because the polar wind supply of mass to large, equatorial distances is 
increased (for fixed M) for large So- Also, as —> 1 (critical rotation) we find 

^(1,0:) = 3(1 -V (21) 

which is > 1 at X > 3^/^/(3^/^ - 1) 3.26. In Figs. 1 and 2, we show ^ and *(5'o, x) versus 
X for various So- Fig. 1 shows the disk density to peak in the range x 1.3 — 2.3 for all So, 
then to decrease rapidly with x for all Sg- 



3. Effect of Gravity Darkening on Disk Extent, Mass, and Emission Measure 

We have seen that poi^) decreases and moves its maximum somewhat to larger x values 
as So increases. However, we need also to consider the extent of the disk, i.e., the lower and 
upper boundaries of equation (17) in x as limited by the magnetic field strength. While 
enhancement of po locally enhances the local contribution per unit volume (~ pj^) to the 
disk emission measure, it makes the material there harder to torque so that the extent of 
the disk is modified. In particular, for example, equation (20) is not valid in practice since 
the rapid decline in B{x), as x goes up, limits the torquing to a finite distance. 

To estimate the effect of including rotational gravity darkening on observable disk prop- 
erties, we need to assess the effect on the inner and outer disk boundaries. Here we do 
so using a somewhat simpler treatment than that in MTD, namely what MTD termed the 
'switch approximation'. In this, the disk is taken to be rigidly torqued by the magnetic 
field (i.e., v = VgX, where Vg = So\/GM/R) out to the distance where the magnetic energy 
density ^ falls below the rotational kinetic energy density Uke = \pdv'^- We have then, by 
equation (9), 

and by equation (17) 

Uke-^Po^SoX (1--) 1 _ 2^2/3 ■ (23) 

So the outer disk boundary x — XouteriSo,^) is given by setting Ub — Uke- Thus Xouter is 
the solution x to 

x' 1-2SJ/3 Sr ^ ' 



where 
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is a measure of field energy compared to the disk gravitational energy. The inner disk 
boundary in the present approximation is simply the Keplerian rotation distance (c.f., MTD) 

•^mner — >~'o • v / 



4. Results and Discussion 

In Figs. 3 and 4, we show Xinner{So) and Xouter{So) versus So for various 7 values, and 
the corresponding colatitudes, on the stellar surface, of Xinner and Xguter in terms of equation 
(13). It turns out that these boundaries do not change greatly with So once So is larger than 
0.2 — 0.3, but change a lot with 7, and that the mass flux reaching the disk comes from a 
rather small range of colatitudes (e.g. for So — 0.6, 45° ^ ^ 70" with 7 = 6 — see Fig. 4). 
Mass flow from the pole (small 9) leaves the star as part of the wind while equatorial flow 
(large 9) does not achieve Keplarian speed. 

According to equation (25), 7 is determined by the magnetic and gravitational fields. In 
order to make mass flux channeling and torquing possible, 7 has to be substantially greater 
than unity. In terms of observations, the magnetic fields of Be stars are no larger than 
hundreds of Gauss. Hence, 7 should probably be in the range of 1 — >■ 10 for Be stars to meet 
this requirement. In Figs. 3 and 4, we also see that for smaller 7, a larger So is necessary 
for a disk (clearly, the outer radius of the disk must be larger than the inner radius). These 
two figures also show that the gravity darkening has a small effect on the outer radius of the 
disk, which is within ~ 5R for appropriate 7 and So- Similar outer boundaries have been 
derived for some stars using different disk models by Dougherty et al. (1994) and Cote et al. 
(1996). 

The detection of disks by polarization, infrared emission, and emission line strength is 
related to their mass and their emission measure which are proportional to Jy pr>dV and 
fyPodV, respectively, where V is the disk volume. If the disk has thickness H{x) — h{x)R 
at distance x then it contains a total number of particles 

271 /"Pouter (So, 7) 

N = / pD(x)h(x)xdx, (27) 

^ ''-'inner \'~'o J 

and has emission measure 

EM ^ —l pUx)hix)xdx, (28) 

where m is the mean mass per particle. Using equation (17) these yield 

^^^'^ 1^^^' - P'^^ - ^°t' - - In'^'dx, (29) 

*^ dinner {,^0) ol 



-8- 



EM = EM^ / TT^JtWiI - - Soil - (1 - IriVdx, (30) 



where A^o = 27ri?^^, and 

'■xo.te.(So,7) x-'^h{x) 1^2/3r-, o2r-, ^ 

n(So) 

where £;M<, = ^-nR^i^f. 

Following the Brown & McLean (1977) formulation, the scattering polarization is P = 
r(l — 3r)sin^i, where r is optical depth, V is the shape factor of the disk and i is the 
inclination angle. Assuming the disk to be a slab with constant thickness H = Rh and 
including the finite source depolarization factor D = a/1 — R^ /r"^ = ^/l — (Cassinelli 
et al. 1987; Brown et al. 1989), then we have the optical depth r. 



T^-^ / n{r, ij)D{r)dr djj, 



16 



PXonter rfl ^ 

/ / n(x,z)D(x)— -dxdz 

Jx Jo ^ + Z 

inner " 



^To x%l--f j-—f- W 1 - — arctan -o?a;, 31 

L. a; 1 — 26^/3 \ x^ x 

^ ''-■inner C / 

where Thomson cross section, n = Pd/i^ is the electron density of the 

disk, and /i is the cosine of the angles between the incident light to the disk and the rotational 
axis. As in MTD, we neglect the absorption and suppose a fully ionized disk. The shape 
factor r yields 

J;^|:^n{r,^,)D{r)^,'drd^^ 



j;^i:;n{r,f^)Dir)drd^^ 
rxonter j'^;') ('j;') arctan -dx 



(32) 



Substituting equations (31) and (32) in the original polarization expression and after some 
reduction, yields the polarization, P, with gravity darkening effects, 

P = Polp (33) 
where Pq — To— ^"1^° and Ip is the integral 



^outer 



Ip= j ^-\l-- r i_2^V3 (34) 

Then Ip — P/Po is found numerically and shown in Figs. 7 and 11 for /i = 0.5 with respect 
to typical half opening angles of the disk about 10", (e.g., Hanuschik (1996); Porter (1996)), 
and the inchnation angle i — 90°. 
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In Figs. 5 and 6, we show the results of equations (29) and (30) for N{So,^) and 
EM(So, 7) together with those obtained when gravity darkening is ignored (using same 
switch approximation). For the latter, we use the integrands as in equations (29) and (30), 
but with So = 0; the same lower limit Xinner = So"^^^ as given by equation (26); and the outer 
limit the solution to equation (24), with 6*0 = on the left. We see that increasing So from 
zero results in a rising disk mass and emission measure up to a broad maximum at 5*0 < 0.5 
and falling back almost to as (So —> 1. It is not surprising to see that gravity darkening 
has strong effects on total number of particles and emission measure, since gravity darkening 
effects significantly reduce the mass flow from equatorial stellar regions into the disk. We 
also note that N/No and EM/EMg ratios have peaks at about So ~ 0.5 for all 7, while for 
no gravity darkening they essentially keep increasing. Similar results are shown in Fig. 7 
for the polarization which depends mainly on the disk electron scattering mass. The more 
the total number of particles in the disk, the stronger the polarization. If gravity darkening 
were neglected, one would get a large Ip so a small Pq, for a given observed polarization 
value P (equation 33). This would imply a smaller po since Po oc po, which implies an 
underestimation of the mass loss rate M (equation 18), if we ignore gravity darkening. 

The previous treatment is for fixed (3=1. In order to see the influence of the velocity 
law on the results, we tried various /? values for a given 7. Fig. 8 shows that wind velocity 
law has minor effects on the disk boundaries and so does gravity darkening. Figs. 9, 10 and 
11 show that slower winds (i.e., bigger /3 values) will lead to much smaller total number 
of particles, emission measure and polarization of the disk, as well as gravity darkening 
significantly decreases these disk properties. From the plots we see that, for small P = 0.5, 
the total number of particles, emission measure and polarization peaks shift slightly to larger 
rotation rate {So ~ 0.6), while a statistical study of observation data indicates the most 
common So ~ 0.7 (Porter 1996) which, in our interpretation, would favor small f3, i.e., fairly 
fast acceleration of winds from the stellar surface. If individual disk dctcctability peaks for 
So ~ 0.6 and actual detection peak for So ~ 0.7, either there is a bias/selection effect in 
operation (Townsend et al. 2003), or there in an upward trend in the frequency distribution 
of So values. 

5. Comparison of MTD with MHD Simulations and Other Work on Magnetic 

Channeling 

It was noted in Section 1 that the MTD model is phenomenological and parametric, and 
not a full solution to the physics equations. It is aimed, like all such models, at describing the 
main features of a system accurately enough to reproduce the essential physics but simple 
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enough to facilitate ready incorporation of additional effects (such as gravity darkening) and 
comparisons with data. It is of course important to evaluate how well the MTD model 
describes the reality when compared with more complete solutions. A full and detailed 
comparison is beyond the scope of this paper but we summarise here the present status as 
we see it of the relation of MTD to recent analytic and numerical work on closely related 
problems. 

One of the earliest studies of the problem which found disk formation was that in a 
'magnetospheric' context was by Havnes & Goertz (1984). Keppens & Goedbloed (1999, 
2000) carried out numerical simulations of magnetised stellar winds with rather weak fields 
and found disk 'stagnation zones' in the equatorial plane. Matt et al. (2000) studied non- 
rotating winds in dipole fields and found persistent equatorial disk structures around AGB 
stars though with a steady throughput of mass leaking through the disk. Maheswaran (2003) 
conducted a detailed analytic study of the MTD situation and found that persistent disks are 
formed for quite small fields though he obtained somewhat tighter constraints than MTD on 
the relevant regimes of magnetic field and spin rate. In the MHD simulations of isothermal 
flow driven outward from a non-rotating star with dipole magnetic flelds, ud-Doula & Owocki 
(2002) found that the effect of magnetic flelds in channeling stellar winds depends on the 
overall ratio of magnetic to flow kinetic energy density, (as did MTD) and obtained disk 
results with rather low fields. 

In contrast to all of these, in the ud-Doula & Owocki (2003) conference paper, based 
on the same code as ud-Doula & Owocki (2002), the interpretation shifts somewhat and 
seems more negative about disk persistence. Owocki & ud-Doula (2003) added rotation to 
their earlier work and concluded that no stable disk could form, with matter either falling 
back or bursting out after a modest number of fiow times. In fact, the MHD code they used 
is incapable of handling the larger fields which MTD argued were required and which are 
recently in fact observed (several hundred Gauss) in some Be stars, so their numerical results 
arc not that relevant. Furthermore, given that observed fields are strong enough so that the 
wind is bound to be steered and torqued to the equatorial regions, if the behaviour of that 
matter were highly unstable as Ud-Doula and Owocki's simulation results, we should observe 
very frequent Be star disk disruption, which we do not. In fact we see no reason why material 
should fall back, given that it is centrifugally supported, until the dipole structure fills up. 
This takes a very large number of fiow times (many years) - hence the the quasi-steady 
formulation in MTD. In the case of weak fields and low rotation, fall back and burst out of 
matter is not altogether surprising, but it is not clear why the Owocki Sz ud-Doula (2003) 
simulation results conflict with those of others. For stronger flelds using order of magnitude 
scaling estimates, they found that disks, essentially like MTD can form and have gone on to 
develop scenarios for strong flelds closely akin to MTD under the names Magnetically Rigid 
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Disk and Magnetically Confined Wind Shocked Disk. They were, however, dismissive of the 
relevance of this to Be-stars. This was not on the grounds of the physics of MTD but over 
the issue of whether a semi-rigid disk near co-rotation can be reconciled with observations 
of Be-stars, specifically spectrum line shapes and the long term V/R variations. The work 
of Teller et al. (2003) suggests that the former is not a serious problem. The issue of the 
V/R variations was emphasised in the original MTD paper which recognised that, if MTD 
is the correct description of Be disk formation, a close examination is required of how the 
V/R variations could arise. At first sight it would seem that the conventional interpretation 
in terms of spiral density waves (induced by the non-spherical potential) in a Keplerian disk 
would not work if the field controlled the disk and another interpretation would have to be 
found but it depends on the rate of viscous diffusive redistribution of disk matter toward 
Keplerian (Maheswaran 2003). Until further testing is carried out we are therefore of the 
view that the MTD model remains a good basic scheme for further modelling work. Wc 
also note that, when MTD is applied across the range of hot star spectral types, it offers a 
remarkably good explanation of the narrow spectral range where disks are in fact detected. 
No other model offers any explanation of this observation. 



6. Conclusions 



We have discussed the phenomenological magnetically torqued disk (MTD) model of 
Cassinelli et al. (2002) for hot (particularly Be) star disk formation in relation to other 
work on magnetically steered wind creation of disk like structures, concluding that, for 
moderate fields comparable to those observed, the description is physically realistic but that 
further work is needed to see if its disk velocity structure can be reconciled with observations 
including hne profiles and V/R variations. We have recognised that the basic model did 
not recognise the effect of spin induced gravity darkening on the latitudinal distribution of 
wind flow and consequently on disk density structure. We have included this effect and 
found that, although increasing So from zero favors disk formation, at high So the polar 
shift in mass flux results in decreasing disk detectability by emission or polarization. The 
fact that detectability (say above half of the height of the peak values of emission and 
polarization in Figs. 6, 7, 10, and 11) covers quite broad rotation range in So, namely 
So ~ 0.25 — 0.80, is generally in good agreement with the fact that Be star rotation rates 
are typically estimated to occur most frequently near 0.7 (e.g.. Porter (1996), and references 
therein; Yudin (2001)), which tends to favour fast acceleration velocity laws. Overall this 
means that in the MTD model, the most easily observable disks are, contrary to naive 
expectation, not expected from the fastest rotators, but from moderate ones, as observed, 
though we note the Townsend et al. (2003) comments on the effect of gravity darkening in 
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inferring underestimated line width rotation rates. Clearly the whole MTD scenario needs 
further work to test it thoroughly, including reconciliation of the phenomenology, numerical 
MHD, and analytic MHD theoretical treatments and work on further diagnostics such as 
X-ray emission from the MTD deceleration shocks of Be stars for comparison with ROSAT 
and other X-ray datasets. 
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Fig. 1. — Disk density ^ vs equatorial distance x for various spin rates, So- Disk density 
increases rapidly then decreases sharply again. The peak value decreases and moves out as 
So increases, but the peak is in the range x ~ 1.3 — 2.3 for all So- 
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Fig. 2. — Ratio of disk density to that without gravity darkening \1/ = ^^(q°^) vs equatorial 
distance x for various spin rates So- 




Fig. 3. — Disk boundaries Xinner and Xouter vs So for various 7. The solid line corresponds to 
the inner radius Xinner, and the other lines to the outer radius of the disk for various 7. Each 
pair of curves with the same line type gives results with gravity darkening (upper curve) and 
without gravity darkening (lower curve). For larger 7, one requires smaller So in order for 
a disk to form {xouter > Xinner)- This implies that for small magnetic fields, a relatively fast 
spin rate So is necessary, while for large magnetic fields, moderate So is adequate. Comparing 
curves with and without gravity darkening shows that gravity darkening has only a minor 
effect on the outer radius of the disk for any given 7, -So. 
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Fig. 4. — Stellar surface colatitudes 9° corresponding to the inner and outer disk boundaries, 
vs So for various 7 . The solid line corresponds to the inner boundary, and the other lines to 
outer boundaries of the disk. Each pair of curves with the same line type gives results with 
gravity darkening (lower curve) and without gravity darkening (upper curve). The possible 
range of colatitudes channeling matter are angles between the curve of solid line type and 
the other curve for a given 7. 
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Fig. 5. — Ratio of total number of particles in the disk relative to arbitrary reference 
value No, vs So for various 7. Lower and upper curves with the same line type are for gravity 
darkening and no gravity darkening, respectively. Gravity darkening significantly decreases 
the results as So becomes large, and produces a peak. For larger 7, a smaller So gives rise 
to the formation of a disk, so the curve becomes wider. The larger 7, the bigger 




Fig. 6. — Ratio of emission measure relative to arbitrary reference value EMo, vs So for 
various 7. Lower and upper curves with the same line type are for gravity darkening and no 
gravity darkening, respectively 
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Fig. 7. — Ratio of polarization -p- relative to arbitrary reference value Po, vs So for various 7. 
The inclination angle is assumed to be 90" (edge-on observation). Lower and upper curves 
with the same line type are for gravity darkening and no gravity darkening, respectively. 
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Fig. 8. — Disk boundaries Xinner and Xouter vs So for various /? with 7 = 6. The sohd hue 
corresponds to the inner boundary Xinner, and the other hues to outer boundaries of the 
disk. Lower and upper curves of the same hne type are for gravity darkening and no gravity 
darkening, respectively. 
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Fig. 9. — Ratio of total number of particles jj- relative to arbitrary reference value No, vs 
So for various (3 with 7 = 6. Lower and upper curves of the same line type are for gravity 
darkening and no gravity darkening, respectively. 




Fig. 10. — Ratio of emission measure relative to arbitrary reference value EMg, vs So 
for various (3 with 7 = 6. Lower and upper curves of the same line type are for gravity 
darkening and no gravity darkening, respectively. 
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Fig. 11. — Ratio of polarization relative to arbitrary reference value P^, vs So for various 
(3 with 7 = 6. The inclination angle is assumed to be 90" (edge-on observation). Lower 
and upper curves of the same line type are for gravity darkening and no gravity darkening, 
respectively. 



